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TECHNICAL MEMORANDUM NO. 387. 

APPROXIMATE CALCULATION OF THE 
STATIC LONGITUDINAL STABILITY OF AIRPLANES.*' 
By The odor Bienen. 

By the static stability of an airplane is- meant its ability 
to withstand forces and moments, which tend to disturb its state 
of equilibrium, by opposing greater forces and moments (whose 
production is the task of the tail group) , and thus to retain 
its original equilibrium. 

In considering static stability, no attention will be paid 
to the motions of the airplane, which are produced by the dis- 
turbing forces, as these come under the head of dynamic stabil- 
ity. Static longitudinal stability has been exhaustively treat- 
ed in numerous articles, the first of which was published in 
Germany in !, Flugsport , 11 1910, by.H- Reissner.* 

The accurate calculation of static stability, especially 
for multiplanes, is now very troublesome- This is especially 
noticeable, when such details as stagger and decalage (inter- 
inclination) are taken into account in calculating the mutual 
effect of the wings. It is a question as to whether such a de- 

* "Eine ^infache Metftode zur angena-iexten Berechnung der stat- 
isc'nen Langs stab ilitat von Ei.n- "and I'oppeldeckern, " from Zeit- 
schrift fur Flugtechnik und Motcrluf tschif f ahrt , July 28, 1926, 
pp. 299-305- This treatise was taken from the course in aviation 
at the Aachen Technical High School, which explains its didactic 
character. 

** See Fuchs-Hopf, "Aerodynaaik," Berlin, 1922, p. 310 ff., and the 
bibliography on p. 459 of the same book* 
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tailed calculation is appropriate when ether influences of fully 
as much importance, are more or less negledted. Such details 
include* 

Effect of propeller slip stream on wings and tail group; 

Deviations of the profile or wing section coefficients in 
experiments with models and full-size airplanes, resulting from 
different Reynolds Numbers and from differences always existing 
between the model and the full-size wing; 

Differences between the assumed and actual lift distribu- 
tion, not only between the upper and lower wings, but also 
throughout the span; 

Differences between the computed and actual position of the 
center of gravity; 

Inaccurate assumptions on the magnitude of the structural - 
resistance, and the point of application of its several coeffici- 
ents, especially at different angles of attack* 

It seems at least desirable to have some simple method for 
calculating quickly and with sufficient accuracy • 

1. The correct position of the center of gravity; 

2- The requisite tail-group dimensions; 

3. The course of the wing and tail-group moments. 

In our deductions, we will first replace the biplane (dis- 
regarding the effect of stagger, decalage and induced drag) by 
an equivalent monoplane, whose dimensions and position in space 
can be approximately determined in a simple manner. 
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In the computation of the balancing of the moments and in 
accord with the suggestion of Von Karman, at whose request we 
employed the following method, we adapted, as the point of ref- 
erence for the moments, the intersection of the wing, chord with 
the projection of the leading edge on the chord, i.e., the point 
to which the moment coefficients of experiments with models are 
generally referred. This resulted in certain simplifications, 
since we did not have to convert the lift and drag coefficients 
into normal force and tangential-force coefficients, in the 
well-known tedious manner, but could utilize directly the c m 
values of the experiment with a model. 

Replacement of the Biplane "by a Monoplane 

We start with a biplane (with stagger and decalage) accord- 
ing to Fig. l and seek an equivalent monoplane whose moment is 

M - CfcE "tE F ges Q (l) 

wherein the moment coefficient c m jp refers to the point of in- 
tersection of the wing chord of the desired monoplane with the 
projection of the leading edge on the chord. We therefore seek 
the momentary effective angles of attack x 0 and x U5 y 0 and y u 
(Fig. l) and the wing chord of an equivalent monoplane. 

For an infinite aspect ratio, we can, in general, write 

where c a t( a }' is the derivative of c a with respect to a y . 
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and a + 6 denotes the direction of the air stream relative to 
the c a = 0 line (Fig.- l) . For normal angles of attack, c a ! 
is a constant, to which we will subsequently return. Hence 

c aS F ges <J = (c ao F o + c au F u) Q 

c a ] (a + 6 - a) F Q + c a ! (a + 6 + a) F u 
°a (a) a " F ges 

from which the effective angle of attack of an equivalent mono- 
plane (at first for an infinite aspect ratio) is found to he 

a = a + 6 - a Jq^S ( 2 ) 
^o + *u 

With a positive deca'lage of a 5 the wing chord of an equiv- 
alent monoplane is inclined toward the mean chord (Mittellage) at 
F - F 

an angle of o . The angle of attack a + 6 is then di- 

* ges 

minished by this araoant, but. with a negative decalage, it is in- 
creased (provided F 0 > F u ) . 

The conversion of the measurement to the correc t- aspect ratio 
is then made in the usual manner. With the most favorable lift 
distribution, we obtain. the angle of attack 

a = oy + A a 



, a = !k x.- f - ) 



The index M denotes the values for the model measurement. 
/Fx 

Hence { rs" denotes the aspect ratio of the model; S F the 
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total wing a.rea of the biplane; b 1 _. the span; while k is a 
function of the ratio, gap to span, and of the upper span to the 
lower span, and of the lift distribution.* 
Thus we ultimately obtain 



= f (a - a 



ges 



2uA 
) 



likewise we can put 



a 



ges 



J 



F - F 1 

The quantity a — ~- ~, which represents the effect of 

* ges 

the decalage, is generally negligible. The other quantities , by 
disregarding the decalage, are reduced very simply to 



t-p 



get 



x Q = x 



y 0 = y 



ges 



sees 



> 



(3) 



J 



If, for comparison, we compute the moment of a biplane and 
of an equivalent monoplane, both referred to the projection of 

the leading edge of the monoplane on the wing chord, we obtain, 

*'Prandtl, "Ergebnisse der Aerodvnamischen Versuchsanstalt zu 
Gottingen," Part II, 1923, p. 15 ff and p. 37 ff. The bending ef- 
fect, which is conspicuous with a relatively small gap, can be 
disregarded in this approximate coinputabion. In like manner, 'we 
can determine the change in the induced drag according to the 
well-known formula r 2 f ^ /tp n i 

and therefrom _ 

c w ~ c wM + ^c w 
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"by omitting unimportant members, 

F u 

t 0 F 0 + c mu tu F u - c ao cosa F Q x 

v x F o 
+ o au cos a i? u = 

r ges 

= OqE (t 0 F Q + t u F u ) . 

We can now write 

c m = c rn ' (c a ) c a + c rnpj 

°a = c a'(a) (« + «)• 

By substituting the corresponding values right and left, 
we obtain 

c m « c a T [ (a + 5 - a) t Q F Q + (a + 5 + a)t u F u ] + 

F Fx 

+ Cmp C*o ? o + % F u3 + Ca' cos a ° Fg g S 2 ° 

F — F 

- [<V c a' (a + 6 - a; ° Fges U ) + ] 

It is obvious that, for o = o, the right and left sides 
agree. If o * 0, we would then have to add to the right side 
(disregarding the higher order of magnitude of the small members) 

v 



2 c a ! o x — * cos a. 

so that the expression for the moment would assume the following 
form: F F 

* E E ^ \ c mE + ^ c a a ' C E P- 2 / 
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The second term in the brackets, however, can obviously be 
omitted, if both the stagger and the decalage do not simultaneous 
ly have large values.- This term is positive when both stagger 
and decalage are simultaneously either positive or negative and 
it then increases the moment- When the signs are unlike, how- 
ever, the expression is negative. 

Calculation of the Static Stability 

We will now turn to the calculation of the static stability 
of a monoplane, to which, as we demonstrated in the preceding 
section, the biplane can be reduced. 

We will first consider the stability in horizontal flight. 
Fig. 2 represents the airplane under consideration- The origin 
of the coordinates is located at the reference point of the mo- 
ments and all quantities are calculated as positive downward and 
backward. The moments are balanced (a = o) at the angle of 
attack i. The angle a refers to the mean-chord line (Mittel- 
linie) and therefore has a different meaning from that given it 
above, where a is the angle of attack of the wing chord. 

S = propeller thrust turning about the center of moments on 
the lever arm s. (3 is assumed to be parallel to the mean- 
chord line); 

• G = total weight of airplane; 

r = distance back to center of gravity; 

h = distance down to center of gravity; 
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W = total structural drag acting on lever arm w (W x w is 
assumed to be independent of a); 

I = distance of center of lift of tail .group from the center 
of moments; 

f = total area of horizontal tail planes. 

Then the total moment of the airplane, referred to the nor- 
mal center of moments of the wing section, is 

k - q F t | 0 m - C w ge3 t a C a - 

+ °ws ? + T F ( C0S(X ~ ¥ sina )} - < 4 > 
In this formula, = the contribution of the wing, 

Cw ctpq *f — ~ — = "the contribution of the pro- 
w S es t cos a * 

peller thrust (in unaccelerated flight), 

c a r CQSa f h gina = the contribution of the 

total weight, applied at the center of gravity of the airplane, 

c ws "| = contribution of the structural 

drag, which, as already mentioned, is assumed to be independent 
of the angle of attack, 

7 -f • \*y > 

— ~r c^i ( ccs a ~ sina 1= the contribution 
x F m \ 1 

of the horizontal tail plane (c n j f = normal force coefficient 

times the area of the horizontal tail plane). 

We now transform . the expression for the tail plane and com- 
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pute the angle of attack i 1 (Fig* 2) in such manner that the 
tail moment t vanishes for a = 0 (hence for the balanced flight 
condition) . 

At first we can disregard — sin a in comparison with 
cos a . Furthermore, approximately 

• c nl = c al cos a = V(a) (a + V + A a F) cos a - 

c a ' is a constant for the existing angle-of-attack region. 
With 'sufficient accuracy, as demonstrated by the experimental 
results, we ca.n write 

c a ' (a) = - : 'V (5) 

1 + 1 . 67 Jg- 
b^ 

A Oy is the downwash angle produced , by the downwash from the 
wings, at the position of the. horizontal tail plane. It is 



F " TT f Vb 2 / 



1 °a 



The index F shows that the aspect ratio of the wing is meant - 

If the airplane is to be balanced for 6 (.Fig* 2), with a 

* This downwash angle was produced by elliptical lift distribu- 
tion resulting from the descending vortex, while the still less 
important effect of the supporting vortex and also the effect of 
the distortion ( rolling up; of the vortex band^ which occurred 
at large c a values, ms disregarded. See Fuchs-Kopf "Aerodynam- 
ik," Berlin, 1S22, p. 145, and especially, H. B. Helmbold rf l, Ueber 
die Berechnung des Abwindes h inter einem rechteckigen Flugel," 
Z.F.HE. 1925, p. 291. Our formula, aside from the disregard of the 
supporting vortex, -represents an upp^r limit for the downwash 
angle. It is sufficiently accurate for practical purposes. 
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symmetrical cross section of the horizontal tail plone, we ac- 
cordingly have 

After inserting: this, we have 

= c a ! 7 a cos a (l.- c a ! 7p) (7) 

Herein 

n \ - 5. 35 . i - 2 / F > 5*35 / _ \ 

°a I 74 \ > c a * l—, /F \ (5a) 

1 + 1.67 n *^ 2 T 1 + 1.67 lK] 

The indexes i and F refer respectively, to the tail and 
wing. We accordingly have 

M _ n _ s 1 _ n r cos a - h sin a . w 

Ft ^ - § es t cos a ° a t q ^ t 



q F t 



I f 

+ t F & (l - ca 1 ?) cos 2 a (4a) 



Taken in order, the terns on the right side indicate* 

The moment coefficient of the wing, 

" " 11 . u n propeller thrust, 

11 " " produced by the weight of the air- 

plane, 

" 11 n of the structural or parasite drag, 

" M » » " tail. 

The moment coefficients of the propeller thrust and structural 
drag are generally of no importance. 

The airplane is to be balanced for the angle of incidence 
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i (Fig. 2), hence for a = 0. For this cr.se we have 

q F t ^ w S es t °G t Ct7S t 

From this we obtain the expression for the requisite dis- 
tance aft of the center of gravity 

~ Cn;c - c w g es c JL + c ws ^ 

f = 5 *- V (8) 

x °ac 

wherein the values c> :a0 , ^w ges o> .and c ao are to be inserted 
for the corresponding angle of incidence I (a = 0). 

All the quantities -^re now known for calculating the moments, 
excepting if for the tail moment. In order to find this, we 
add the other moment coefficients and plot them against a or 
c a . Then if is so determined that, according to the degree 
of stability desired, the moment coefficient of the horizontal 
tail plane for the existing a- or c a is equal to or greater 
than the maximum wing or. other moment. 

Attention is hereby called to the fact that s and h are 
calculated positively downward from the reference point of the 
moments (Fig. 2) and are therefore partly negative on biplanes 
and low-wing monoplanes.. It is known that the small angle of 
attack combined with large drag coefficients in some airfoils 
(e.g., those with leading edges sharply curved downward) can be 
dangerous for low-winged monoplanes, is there is not a suffici- 
ently large stabilizer provided. 
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Our formula needs to be supplemented for steep gliding 
flight with engine stopped. In this case, the weight of the 
airplane can be divided into one component in the direction of 
the lift and another in the direction of the drag (Fig. 3). In 
unaccelerated flight, these components are respectively equal 
to the lift and drag, but in opposite directions- In equatioh 
(4a) for the total moment, the expression 



ges t cos a 

is then to be replaced by 

r sin a + h cos a . 
c w ges t 7 

we must therefore find whether' 



> £ sin a + £ C os a 



t cos a < t 
The moment formula would then read' 

^ = c m - c, 7 ges -I Blna^ h cos a _ Ca r coscc - h sina 

+ c ws * + | | ca 1 I (1 - ca ! F) cos 2 a (4b) 

This is of no practical importance, since the moment coeffi- 
cient produced by the drag is generally small in comparison with 
the other quantities and, besides, s/coscc can hardly be smaller 
than r sina + h cosa. Moreover, the last equation (4b) also 
applies to engineless airplanes (gliders). 

The expression for the requisite aft position on a glider 
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can be found by putting, as before, a = 0. We then helve 

c m0 - c w ges o t + c ws 2L 

x = r— 1 £ (8a) 

x c ao 

In closing, we wish, to deduce one more simple formula for 
the approximate calculation of the dimensions of the horizontal 
tail plane. For this purpose, we write the equation for the 
equilibration of the moments in the following form: 

= c * ! (ca) 0 a'(a) (a + 6) + %p - ( %LjL + c wS j | - 

- c a T (a) (a + 6) (| + a fj + I f c a ! j a (l - c a ' F ), (9) 

in which we develop sina and cos a and disregard the terms 
with higher powers of a. 

• The stability formula 'now reads 

ff > 0 (10) 

or, expressed in words, if we are to have stability, thenj with 
increasing angle a, a more rapidly increasing positive (hence 
nose-heavy) moment must be produced and, conversely, with a neg- 
atively increasing angle, a, a greater falling negative (hence, 
tail- heavy) moment must be produced. Or, more briefly, a re- 
storing moment must be produced in disturbances of the state of 
equilibrium. 

We differentiate equation (9) with respect to a and obtain 
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_B L = c , [c . _ | _ Hl! s h (2 a +6)] 

3 a q F t a l -m t n t T t x 

If we now disregard the terms containing s/ t and h/t, 
and write approximately for r/t the expression c^qI Cao -'equa- 
tion (8), we finally obtain 

oa'pf 5aft- c m > (ca) ) 

l *> A£ao y 

t c a »l (1 - c a ' F ) 



An Example 

On account of cur many omissions, we considered it advisa- 
ble to test, by an example, the resulting differences, as com- 
pared with the accurate method of computation. For this pur- 
pose, we chose the example published by the D.V.L. (Deutsche 
Yersuchsanstalt fur Luftfahrt). The dimensions of the biplane 
considered are represented in Fig. 4, and are numerically as 
f ollows: 

Upper Wing Lower Wing 



Area F Q =19.2 m 2 (206-7 sq.ft.) F U =18 -0 m 2 ( 193.8 sq.ft.) 

Span b o =!2.0 m ( 39.37 ft.) bu=13.0'm ( 39-37 ft.) 

Mean chord t 0 = 1-6 m ( 5-25 ft.) t u = 1.5 in ( 4.92 ft.) 
Airfoil Gottingen 398 Gottingen 398 

Distvaft- r 0 = 0.82 m( 2-69 ft.) r u = 0.24 m( 0.79 ft.) 

Diet., down h 0 = 0.93 m( 3-05 ft.) h u = 0.67 m( 2.20 ft.) 
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Decalage a = i© 

Stagger p = 20° 

Gap h = 1.6 m (5.25 ft.) 
Elevator- 

Area f = 3.8 m 2 (40.9 sq.ft.) 

Span b H = 3-74 m (12.27 ft.) 

Chord t H = 0.95 m ( 3.12 ft.) 

Distance of center of gravity 

from elevator axis, I = 4.44 m (14.57 ft.) 

Angle of incidence of elevator 

to mean chord of wings , . a = 0.4° 

Ooef . of structural drag > c T;s = 0.02 

Lever arm of structural drag, 0.15 m ( 0.49 ft.) 

Height of propeller axis above 

center of gravity, 0.2 m (0.66 ft.) 

The wings are balanced at an angle of attack of 4.5° and 
c a = 0.77. If the position of the equivalent monoplane is calcu- 
lated according to the above formulas, the distances from the 
leading edge to the center of gravity are 

Aft r - 0.553 rn (l.81 ft.) 

Down h = 0.150 m (0.49 ft.) 
According to equation (8), a distance aft of r = 0.558 m (l.83 
ft.) would be required. The angle of incidence of the elevator 
to the mean chord of the wings is found according to equation 
(6) to be o = 0.6° (against C, 4°) . 

In Fig. 5, the airplane polar is plotted accurately and the 
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points are also given (+) for the approximate polar, disregard- 
ing the stagger and decalage- (The "bending effect" is also 
disregarded here.) 

Fig. 6 shows the course of the moment curve both by the ac- 
curate method and by the .approximate method- The' moment curves 
naturally differ, for one refers to the leading edge of the 
equivalent monoplane and the other to the center of gravity of 
the airplane. In both crses, the wing moment is zero, with al- 
most the same angle of attack. 

On the basis of the preceding considerations and of the ex- 
ample, I think it is not too much to claim that the simple meth- 
od shown is entirely satisfactory for practical purposes. The 
deviations in comparison with the accurate method of computation 
fall within the limits of the errors in computing with a slide 
rule. 

Summary 

A decalaged and staggered biplane can be approximately re- 
placed by a monoplane whose aspect ratio is given by k — §- • 
The moment of this monoplane is 

M = c mE t S F ges q (i) 

in which c^e is the moment coefficient of the given airfoil 
calculated for the above aspect ratio, 
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The airfoil, v;ith positive stagger a toward the mean chord 
(Fig. l), is negatively inclined rith the angle 



J^Z^Jk , (2) 



*c;OS 



while the eft and depth position of the equivalent monoplane 
with the notations of Fig- 1 is given by' the expressions 



= x - ^ 



- ges 
F 



^ (3) 



y 0 = y v^ - • 
° *ges s 



The course of the total moment for an airplane is given by 
the equation 

M - ^ „ e n r cosa - h sin a , n w 

ofTl ~ ^ ' 7 Ses t cosa ° a t ws t 

+ | | c a 'i cc (l - c a ' F ) cos a a , (4b) 

the last term on the right being the coefficient of the tail mo- 
ment, in which 

, _ 5-25 . N 

Pal TT^T (o) 

1 + 1.67 (V ' 



■b ' 



I 



- , 2 / F \ _ o-2 5 , _ , 

n V ° 'F l + 1.67 V 

-b <F 

The aft position of the center of gravity is given by the 
equations 
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c nio - Cy/ ges o t + c^rs ■? 
- - ■ — - (engine driven (3; 



•o airplane) 



and 



w 

Crno - c^r ges if + °WS T , 

f = 1 (Sngineloss (8a) 

x c ao airplane) 

From the* stability condition, we obtain, for the size of 
the elevator, appro:: j mately' 

if. . ^v°a.o (n) 

* cxl c a i j a - c a i F ) 

Supplement 

Professor Von Karman called my attention to the following 
method by which allowance can be made for the effect of the 
stagger on the stability. 

Since the effect is only slight, we will be satisfied with 
an approximate calculation, assuming, as the basis of our deduc- 
tions, that the upper and lower wings do not differ very great- 
ly from one another, either in their dimensions or in the con- 
ditions of flow to which they are subjected. 

We will consider the effect of the lower supporting vortex 
on the upper. The circulation of the lower vortex follows from 
the. expression 

I T u v bu = 7 £ °au F u 

to 



S 2g 



c F 
°au 1 1 

2" " 



p _ °au x u 
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A vertex of the strength T generates, . at the distance r, 
a velocity of the magnitude 

A v = — E — 5 
2 '.i r 5 

perpendicular to the radius vector. 

In our case, the following speed increment is added for 
the upper wing, as the resalu of the circulation about the 
lower Tiring. 

T 1 c F 

L v ou 3 tt r ■ <i it bu r* 

As a result of this speed increment, the upper wing, with 
a positive stagger of (3 degrees (Fig- 7), suffers a change in 
the angle of attack of approximately 

a A Tr ou O ~ a ) 

A a ou - : w : • 

If we disregard the speed increment, we obtain, a,s a result 
of the change in the angle of attack, the following lift incre- 
ment for the upper wing: 

A A ou = Aa cu Ca - l 7 f- F 0 sin (6 - a) 

= £s& Ell 1 ca J 7 Si F 0 sin (0 - a). 
4 11 D u r 27 

For the lower wing we obtain a reduction of 

■ AA tt0 =|f Is 1 c a .7 || F U s m (p - a), 
in the lift* 
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Under the initial suppositions, r/e can now assume that the 
total lift is not changed and that 

A A ou = - AA i0 

or that 

AA = |a_ i 0a . y Ji, sin (3 - a), 

in which t denotes the chord and F the wing area. 

Due to the change in the lift, there is generated, however ^ 
at positive angles of attack, a tail-heavy moment of 

A M = - A A r sin O - a) 

= - S t c a ' y ~- F sirf (P - a), 
4tt 2 g 

the corresponding moment coefficient '.being 



A c n = = - |- a - c a ' sin^ (3 - a) 

' ! 7 t E F ges 



The angle p. - a may be designated as the effective angle 
of stagger. When the a.ngie of attack is equal to the angle of 
stagger, the effect of the stagger vanishes. It also vanishes 
when the lift becomes zero. The coefficient reaches its maximum 
value at 

a ~ 1 (p - 2 6) , 

in which 6 is again the angle between the line c a = 0 and 
the mean-chord line (Fig. 7) . 

In Fig. 8, the course of the supplementary moment coeffici- 
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snt is plotted against the angle of attack a + 6 for the Gottin- 
gen airfoil 426, with staggers of 0°, +20°, -20°. 

The effect is generally very slight. The tail-heavy moment 
coef f icients , occurring with negative stagger (which is very sel^ 
dom met with) at large angles of attack, are small in compar- 
ison with those of a single- wing (monoplane). 

The following are the necessary formulas for the stability 
calculation, given in the more accurate form. First we have 
for the wing moment of the biplane 

11 = 7 ~t E F ges [CmS + 2 c a < a ~ %^ 2 cosS 



T 



" sin2 0 - a)] (12) 

The second term in the brackets comes from the stagger and 
decalage. Since we can put cos a as approximately 1, this 
term produces only a parallel shifting of the % values. The 
third term is simply a result of the stagger. 

In our moment formula (4a), we would therefore have to 
substitute the bracketed expression from equation (l2) in place 
of c m . 

For the more accurate location of the center of gravity, we 
write f3 



in which a denotes the angle between the wiift direction and the 
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mean chord (Figs-. 2 and 7); 5, the angle between the mean chord 
and the line c a = 0; 3, the angle of stagger (with reference 
to a perpendicular to the mean chord)- Corresponding to equation 
(8) x we then have, for the location of the center of gravity, 

F F c s s 

r c mo + 2 c a' * -%~ \ ~ ~%~ 6 sin2 ? ~ c w ges t + c W s t 



_ = . . .. "S es 

* c ao 



(13) 



Translation by Dwight M. Miner, 
National Advisory Committee 
for Aeronautics.. 
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Figs. 1,2 & 3. 



c a = 0,line N 
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Fig.l Biplane replaced by equivalent monoplane 
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<vi l-ean chord line 
Fig. 2 Airplane measurements and quantities 
required for computation. 
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Fig. 3. Airplone in steep glide. 
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Figs. 4 £ 5. 




Fig. 4 Diagram of wing dimensions. 
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Fig. 5 Course of biplane polar by both 

accurate and approximate computation. 
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Figs. 6, 7 & 8. 
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Fig. 6 Koment curves. 
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Fig. 7 Effect of stagger cn the effective angle 
of attack. 
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Fig. 3 Course of supplementary noment coefficient 
produced by various angles of stagger. 



